AD  A  I  0  4>  \  <\  ? 


AD 


TECHNICAL  REPORT  ARBRL-TR-02361 


A  LINK  BETWEEN  SHAPED  CHARGE 
PERFORMANCE  AND  DESIGN,.. 


I 


Ralph  d Shear 
Frederick  S./Brundick 
John  T.  Garrison 


■  September  1981 


OCT  2  7 


A 


US  ARMY  ARMAMENT  RESEARCH  AND  DEVELOPMENT  COMMAND 

BALLISTIC  RESEARCH  LABORATORY 

ABERDEEN  PROVING  GPOUND.  MARYLAND 


1 


Approved  for  public  release;  dlstributior.  unlimited. 


81  10  27  2?7 


I 


'  •/ 


Destroy  this  report  when  it  is  no  longer  needed. 
Do  not  return  it  to  the  originator. 


Secondary  distribution  of  this  report  by  originating 
or  sponsoring  activity  is  prohibited. 

Additional  copies  of  this  report  may  be  obtained 
from  the  National  Technical  Information  Service, 

U.S.  Department  of  Commerce,  Springfield,  Virginia 

221S1. 


The  findings  in  this  report  are  not  to  be  construed  as 
an  official  Department  of  the  Army  position,  unless 
so  designated  by  other  authorized  documents. 

Th«  ua«  of  ti -ad*  noma  or  manufacture  •  ncmaa  in  this  raport 
dots  not  constitute  indoraamant  of  any  oomareial  product. 


1 1 


UNCLASSIFIED 


SECURITY  CL  A*tl*lCATIOM  O'  TmI*  MOI  TUMI  Dw  BnWQ 


RCFORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1  NKClPKHT'S  CATALOG  NUMBER 

*  TITLE  kMII>j 

A  LINN  BETWEEN  SHAPED  CHARGE  PERFORMANCE 

AND  OESIGN 

TVPI  OF  REPORT  ft  PERIOD  COVERED 

Final 

ft  RCRPORMINO  ORO.  REPORT  NUMBER 

T  AuThoATM 

Ralph  E.  Shear 

Frederick  S.  Brundlck 

John  T.  Harrison 

a.  contract  or  gramt  Humbert*; 

US  Army  Ballistic  Research  Laboratory 

ATTN:  DRQAR-BLV 

Aberdeen  Proving  Ground,  MO  21005 

HHHi 

tt  CONTROL  LIMA  OrflCE  HAMt  AMO  AODAEtt 

US  Army  Armament  Research  a  Development  Command 

US  Army  Ballistic  Research  Laboratory 

ATTN:  DROAR-BL 

It.  RtRORT  OATS 

SEPTEMBER  1981 

it  mumosr  or  r aoei 

27 

i*  monitBWmg  /JttWtlBVI  uMwl  Mmu  Ami  ChwIMn*  OWmJ 

It.  StCURlTV  CLAM  (ml  MM*  NRMtj  ““ 

UNCLASSIFIED 

II  OHriUftuTlO*  tt  At tMfNt  fmi  iMa  Hm»mH) 

Approved  for  public  release;  distribution  unlimited. 

IT  OIITailuTIO*  ITATtMENT  (ml  rnm  aMrail  wNfi  k  Nut  M.  II  rnifmrmml  hmm  Aapw*; 

l»  *EY  BOAOt  (Cmmtlmmm  mm  tmmmtmm  m!4m  II  fimmmmmmrr  mat  imrnilly  *r  M«c*  iwmMkJ 

Shaped  Charge 

Jet  Break-up  Time 

Jet  Parameters 

M  *MT*A€T  itMtM  m  him,  •**  If  mmmmmmtm  «M  irnmmrtlr  *r  Mm*  mMJ 

It  Is  Illustrated  that  the  penetration  performance  of  a  shaped  charge 
determines  *best*  values  of  parameters  In  the  DIPerslo,  Simon,  and  Merendlno 
theory  of  penetration  by  shaped-charge  jets  and  that  It  Is  possible  to  relate 
these  penetration  parameters  to  design  parameters  such  as  cone  angle  and 
liner  thickness. 

DO  .'2Tr»  W1  twnoMor  >novM'SOMolctc  UNCLASSIFIED 


SECURITY  CL  AMiriCATto*  O'  Thi*  fSt  7*Cwa  D«i  Tmwmi 


I. 

II. 

III. 

IV. 

V. 

VI. 


Table 

1 

2 


Figure 

1 


2 

3 

4 

5 

6 
7 


TABLE  OF  CONTENTS 

INTRODUCTION 

DETERMINATION  OF  UMJN  AND  t 
AN  EXAMPLE 

VIRTUAL  ORIGIN  APPROXIMATION 
DETERMINATION  OF  THE  ENERGY  CONSTANT 
SUMMARY 

ACKNOWLEDGEMENTS 

LIST  OF  TABLES 

Some  Calculated  and  Observed  Jet  Data  for  the  3.81  cm  11 

Copper  Liner  Shaped  Charge 

Experimental  and  BASC>Code  Generated  Values  of  Jet  Tip  12 

Velocity  for  Selected  3.81  cm  Aluminum  and  Copper  Liners 

LIST  OF  FIGURES 


Calculated  and  Experimental  Penetration  Depth  vs  Standoff  9 
Distance  from  Virtual  Origin  for  the  BRL  Precision  Shaped 
Charge 

Penetration  Depth  vs  Virtual  Standoff  Distance  for  the  14 

3.81  cm,  Cu,  Conical  Liner,  Cone  Angle  =  20° 

Penetration  Depth  vs  Virtual  Standoff  Distance  for  the  15 

3.81  cm,  Cu,  Conical  Liner,  Cone  Angle  =  40° 

Penetration  Depth  vs  Virtual  Standoff  Distance  for  the  15 

3.81  cm,  Cu,  Conical  Liner,  Cone  Angle  =  60° 

Penetration  Depth  vs  Virtual  Standoff  Distance  for  the  17 

3.81  cm,  Cu,  Conical  Liner,  Cone  Angle  =  90° 

Calculated  Virtual  Origin  vs  Cone  Angle  18 

Calculated  UM^n  vs  Cone  Angle  19 

DISTRIBUTION  LIST  23 


Accession  Tor 

NTIS  GRA4I 

n 

DTIC  TAB 

□ 

Unennounced 

□ 

Justif icfttion_ 

- Pag 

P- 

K  ’  ~*.r4  H:ticn/ 

/.v.v.i-  I'.ir.v 

Co-’.ea 

j  r  i 

/v 

:/or 

20 

21 


3 


w 


I .  INTRODUCTION 


DiPersio,  Simon,  and  Merendino-1  presented  equations  to  determine 
the  penetration  depth  and  hole  volume  associated  with  a  shaped  charge  jet 
impacting  a  given  target.  In  particular,  given  jet  and  target  material 
densities,  pj  and  pt,  jet  break-up  time,  t^,  initial  jet  tip  velocity, 

o  .  -  '  ' '  the  penetration  depth,  as  a 

In 


V. 

J 


1  -  -V 

minimum  penetration  velocity*. 


min* 


function  of  the  virtual  stand-off  distance,  ZQ,  can  be  computed. 


addition. 


if  given  average  jet  diameter. 


and  an  energy  constant,  C, 


DiPersio,  et  al  provide  equations  which  enables  one  to  calculate  the 
hole  volume  associated  with  the  penetrating  jet.  DiPersio,  et  al 

obtained  values  of  t,,  U  .  ,  and  V.°  from  experimental  measurements  for 

l  min  i 

J  Q 

a  precision  shaped  charge,  with  a  42  conical  liner,  and  calculated  the 
total  penetration  depth  as  a  function  of  stand-off  distance  for  this 
particular  charge.  They  obtained  favorable  agreement  with  experimental 
measurements  of  penetration  depth  at  various  stand-off  distances  where 
the  stand-off  distance  is  the  distance  from  the  base  of  the  liner  to  the 
target . 


A  question  raised  by  one  of  the  authors  (J.T.H.)  was,  "Under  what 

conditions  does  the  experimental  penetration  depth  -  stand-off  data  and 

hole  volume  -  stand-off  data  determine  or  infer  the  values  of  C,  U  .  , 

min 

and  tj?",  i.e.,  the  parameters  utilized  in  the  DiPersio,  Simon,  and 
Merendino  (DSK)  equations. 


A  partial  answer  to  this  question  was  given  earlier  by  Majerus  and 

Scott^,  who  utilized  a  modified  form  of  the  DSM  equations  and  investigated 

the  round-to-round  variability  of  C  and  U  .  .  Majerus  and  Scott  provided 

min  * 

a  computational  method  of  determining  C  and  U  ^  from  experimental 

penetration  and  hole  volume  -  stand-off  data.  In  their  method,  they 
required,  in  addition  to  target  and  material  properties,  location  of 
virtual  origin,  jet  break-up  time,  tj,  jet  tip  velocity,  jet  diameter, 
etc. 


/?.  DiPersio ,  J.  Simon ,  and  A.  Merendino,  "Penetration  of  Shaped  Charge 

Jets  into  Metallic  Targets,"  BRL  R-1296,  September  1965,  (UNCLASSIFIED) . 

*  (AD  #476717) 

Also  called  an  interaction  parameter ;  see  Reference  2. 

2 

J.  Majerus  and  B.  Scott,  "CUMIN:  A  Computer  Code  for  Determining 
Certain  Jet/Target  Parameters  from  Experimental  Data,"  ARBRL-TR-02129, 
December  1978 ,  (UNCLASSIFIED) .  (AD  #B035331L) 
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In  the  following,  we  show  that  functions  of  the  DSM  parameters,  t^ , 
Ujnin»  and  C,  can  be  determined  from  experimental  penetration  and  hole 

volume- stand-off  data  or,  in  fact,  from  desired  penetration  performance 
data.  These  functions,  together  with  specification  of  V^0  and  jet  diameter 

dj  yield  estimates  of  tj ,  Umin,  and  C.  Since  Vj°and  d^  are  readily 

determined  from  the  BASC3  code  and  only  require  knowledge  of  material 
densities,  some  explosive  properties,  liner  thickness,  c,  and  cone 
angle,  a,  the  methodology  provided  herein  enables  one  to  calculate  these 
DSM  parameters  without  additional  experimentation.  Such  a  procedure  may 
be  useful  in  shaped  charge  design  problems. 


II.  DETERMINATION  OF  UMIN  AND  t 

Letting  x  =  vP  tj ,  and  y  =  Um^n  tj  then  the  total  penetration  of  the 
jet  into  the  target  is  given  by* 


PT  -  Zc  f  {  x/(l*Y)y  }  1/y  .  i  ] 

whenever 

0  i  Zo  ~  (1+Y)y  t(i+Y)y/x]1/Y 


CD 

(2) 


where  y  =  /p  /p . ,  or  by 

^  J 

pt  -  [Cin)x1/(Y+1)zoY/cl+Y)-/(in)yx1/(Y+1)zov/(ln)]/Y-zo  (3) 

whenever 

(l^-YjyfCi+YJy/x]1^  <_  zq  <_  x  C4) 

or 

PT  ■  [x  -  *'y(x+YZo)  ]  /  Y  (5) 

whenever 

x  1  Z0  <  x  (x/y  -  1)  /  y  (6) 


3  <7.  Earrison,  "Improved  Analytical  Shaped  Charge  Code:  BASC", 
ABBRL-TH-  02200 ,  March  1981.  (AD  #A100275) 

* 

'■'quctionn  (?.?)-(?.$)  of  reference  2. 
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Equations  Cl)  -  (6)  enable  one  to  calculate  the  total  jet  penetra¬ 
tion  as  a  function  of  stand-off  from  the  virtual  origin,  ZQ,  whenever 

x  and  y  are  known.  We  note  from  (2),  (4),  and  (6)  that  the  boundary  of 
each  region  is  also  a  function  of  x  and  y.  Thus  if  x*and  y*  are  known 
values  of  x  and  y,  then  this  specification  determines  a  partition  such 
that  given  a  value  of  ZQ  one  can  determine  the  corresponding  value  of  P^,. 

If  we  are  given  {  (PT  .  ,  Z  .)  }  for  i  =  1,  ....  N  and  where 

1,1  0,1 

P„  .  is  either  the  observed  value  of  PT  at  =  Z  .  or  is  the  desired 
T,i  T  o  o,i 

performance  at  Z  =  Z„  .  then  we  can  obtain  "best"  values  of  x  and  y,  i.e. 

0  0,1 

x*,  y*  as  follows.  We  note  that  the  boundary  between  each  region  of 
validity  for  equations  (1),  (3),  and  (5)  is  a  function  of  x  and  y,  thus 
for  each  value  of  x  and  y,  we  can  compute  the  value  of  P^  =  f  (x,y,ZQ) 

for  any  given  value  of  ZQ.  If  not,  then  the  values  of  x  and  y  lie 

outside  the  feasible  region.  We  let 

N  2 
H(x,y)  =  I  [  PT>.  (Z0>.)  -  f  (x.y,Z0ji)  ]  (7) 

and  we  determine  x*,  y*  such  that 

H  (x*,y*)  <_  H  (x,y)  for  all  x,y.  (8) 

If  Vj°  is  known,  then  tj  and  Umin  follow  from  the  definition  of  x  and  y. 

III.  AN  EXAMPLE 

Experimental  data  for  the  BRL  Standard-Shaped  charge  are  provided 
by  DSM.  Included  within  this  data  are  total  penetration  vs.  stand-off, 
jet  break-up  time,  initial  jet  tip  velocity,  and  minimum  penetration 
velocity  U  .  .  We  have  utilized  the  penetration  stand-off  data  for 
stand-off  m  n  distances  through  20  cone  diameters  (we  did  not  use  the 
penetration  depth  at  25  cone  diameters)  in  equation  (7),  i.e.,  we 
obtained  the  solution  x*,  y*  from  obtaining 


MIN  I  [  pt  ,  (I  j)  -  f  (*.y.z0ii)  1  2  (9) 

x,y  i=i 

from  which  we  found 

x*  =  85.905  cm  (10) 

y*  ■  11.41  cm 
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DSM  reported  that  *  0.830  cm/ysec  thus  since  x* 

and  y*  =  U  .  tf  we  have 
mm  l 


Vj  t? 


tj  =  103.5  ysec 

U*.  =  0.110  cm/ysec 

min 


(U) 


as  compared  to  DSM  experimental  values  of 

tj  =  103  ysec  ^2) 

U  .  =0.10  cm/ysec 

min 

It  is  appropriate  at  this  point  to  recall  that  V^.°  can  be  calculated 

from  the  BASC  code,  thus  the  above  calculation  can  be  performed  without 
knowledge  of  the  experimental  value  of  V^°. 

Since  the  determination  of  x*  and  y*  also  results  in  the  deter¬ 
mination  of  the  corresponding  region  of  penetration,  i.e.,  Z  .  corres- 

0,1 

ponds  to  a  region  in  x  -  y  space,  the  penetration  is  also  calculated  -  and 
required  in  the  minimization  of  (9) .  The  calculated  penetration  vs  virtual 
stand-off  distance  is  shown  in  Figure  1  along  with  the  experimental  values 
of  the  penetration  depth.  The  agreement  is  excellent. 

The  minimization  of  (9)  was  accomplished  by  utilizing  the  "Complex 
Method"  due  to  M.  J.  Box4.  This  method  requires  only  function  evalua¬ 
tions  and  not  derivatives;  thus  the  method  is  ideal  for  this  particular 
application. 


IV.  VIRTUAL  ORIGIN  APPROXIMATION 

In  the  above  example,  the  penetration  depth  was  given  as  a  function 
of  the  virtual  stand-off  distance.  In  the  DSM  report,  the  authors 
obtained  the  location  of  the  virtual  origin  from  flash  radiograph 
measurements;  however,  in  many  other  reports,  the  virtual  origin  is 
either  not  given  or  is  approximated  by  a  "rule  of  thumb".  For  example, 
DiPersio,  Jones,  et  al5  use,  from  past  experience,  the  rule  "the 


'  'x,  'A  V*’”  of  Ton»traine>i  Optimf  nation  aw*  o. 

I'f  cr  Vct'c’s,  '  -'o*puter  J. ,  ^  (1J0S}. 


Comparison 


v.  nivrreio,  7.  Jones,  A.  uervnrf'-no,  <vno  J. 
\-t.  from  i’Z  Cali'  er  f>:  Tpr  *  *  arpev 

iVjw/ve v  ftrrTAS~7FT':V) . 


Simon,  '’fiaraote.rietics  c* 
Sapp#*  ir  '  Aluminum  Liners, 

(AD  #823839) 
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PRECISION  SHAPED  CHARGE 


approximate  location  of  the  virtual  origin  of  a  highly  confined  charge, 

....  is  three-fourths  of  the  liner  height  In  attempting  to 

determine  t  and  U  .  from  the  data  of  reference  5,  we  found  that  the 
min 

above  rule  did  not  result  in  adequate  agreement  between  computed  and 
experimental  values.  Therefore,  we  modified  our  computational  pro¬ 
cedures  and  let 


Zo  =  B  +  S.  (13) 

where  B  is  the  distance  from  the  base  of  the  liner,  along  the  cone 
axis,  to  the  apparent  origin  of  the  jet,  and  S  is  the  stand-off 
distance.  Thus  equation  (7)  becomes 

N  , 

H  (x.y.B)  =  E  (  P  (S  B)  -  f  (x,y,B,S  )  ]z  (14) 

i=l  1 

so  that  we  now  seek  x*,  y*,  and  B*  such  that 


H  (x*,y*,B*)  <  H  (x,y,B)  . 

Utilizing  the  penetration  data  of  reference  5,  equation  (14)  was 
minimized.  In  this  minimization  process,  we  constrained  B  to  lie  in 
the  interval 


where  B  *  height  of  cone  +  distance  allowed  for  liner  retainer  ring 

(:  1.4  cm).  For  the  20°,  60°,  and  90°  conical  liners,  the  resulting 

agreement  of  calculated  and  experimental  jet  break-up  time  was  excellent. 

For  the  40°  copper  liner,  we  found  that  if  B  was  taken  to  be  twice 

max 

the  liner  height,  then  good  agreement  could  also  be  attained  for  this 
case.  In  Table  1,  we  present  calculated  break-up  times  t^*  and  observed 

values  tj,  calculated  minimum  penetration  velocity  and  the  calculated 

location  of  virtual  origin  B*. 

In  obtaining  the  jet  break-up  times,  tj*,  listed  in  Table  1,  we 
used,  in  each  case,  the  corresponding  experimental  jet  tip  velocity 
reported  in  reference  5;  however,  it  is  noted  again  that  the  jet  tip 
velocity  can  be  calculated  from  the  BA.SC  code  utilizing  liner  thickness, 
e,  apex  angle,  a,  and  explosive  and  liner  material  properties.  In 
Table  2,  we  compare  the  BASC-code  generated  values  with  the  experimental 
values  for  some  of  the  3.81  cm  copper  and  aluminum  liners  of  reference  5. 
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Table  1.  Some  Calculated  and  Observed  Jet  Data  for 
the  3.81  cm  Copper  Liner  Shaped  Charge 
(asterisk  denotes  calculated  value) 


'Sj,  c, 

* 

t, ,  nsec 

t, ,  nsec 

U*.  ,  cm/nsec 
mm 

* 

B  ,  cm 

20° 

41.5 

40.8 

0.18 

12.2 

o 

O 

62.5 

63.9 

0.16 

10.0 

60° 

65.3 

66.7 

0.14 

4.4 

kO 

o 

0 

63.4 

64.3 

0.11 

0.0 
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labie  2.  Expei imental  and  BASC-Code  Generated  Values  of  Jet  Tip 
|  Velocity  for  Selected  3.81  cm  Aluminum  and  Copper  Liners 

i 


Cone  Angle 


Material 


\\°,  cm/usrt- 


V  ,  cm/usec  (BASC) 


In  Figures  2-5,  we  have  plotted  the  "best"  penetration  -  virtual 
stand-off  curves  generated  by  minimizing  (14)  for  each  of  the  3.8  cm  (1.5") 
copper  liners  of  reference  5.  In  each  case,  we  used  the  average 
penetration  values  for  each  liner  and  we  have  plotted  these  average  values, 
for  comparison,  on  each  figure.  With  the  exception  of  the  20°  liner, 
the  agreement  is  satisfactory. 

In  Figure  6,  we  have  plotted  the  computed  "best"  value  B*  of  the 
virtual  origin  location  as  a  function  of  cone  angle  for  the  3.81  cm 
copper  liner  and  the  42°  BRL  precision  shaped  charge  of  reference  1. 

The  plot  indicates  that  the  virtual  origin  location  is  approximately 
linear  with  respect  to  cone  angle. 

* 

Finally,  in  Figure  7,  we  have  plotted  "best"  values  of  ^  as 

a  function  of  cone  angle  for  the  1.91  and  3.81  cm  copper  conical  liners 
of  reference  5.  It  appears  that  for  the  liner,  explosive,  and  target 
complex  of  reference  5  that  U*.^  is  approximately  linear  with  respect  to 

cone  angle  and  does  not  depend  greatly  upon  the  cone  base  diameter  for 
these  scaled  liners.  Also,  on  Figure  ~  we  have  plotted  U*.  which 

was  calculated  from  the  penetration  stand-off  data  of  reference  1.  Ke 
note  that  both  the  explosive  and  target  properties  have  changed  for 
this  case. 


V.  DF.TF.RMINATION  OF  TUI.  F.NF.RGY  CONSTANT 


The  hole  volume  produced  by  the  penetrating  jet  can  be  calculated 
for  each  region  of  penetration  by  the  equations  (38),  (40),  and  (42) 
of  the  DSM  report.  For  example,  for  region  1+ 


where 


tt  =  5  x  { 


f 


(WlL-i 

X  1 


F, 


24C 


Pi 


(IS) 


For  each  of  the  other  regions,  each  equation  is  a  function  of  f, ,  x,  ZQ, 
and  y.  We  have  shown  previously  that  x  and  y,  i.e.,  x*,  and  y*  can 
be  obtained  by  minimizing  (7)  or  (14),  and  have  noted  that  dj  and  V. 

can  be  obtained  from  the  BASC-code,  thus  if  we  denote  the  calculated 
hole  volume,  in  its  appropriate  stand-off  region  by  g  (Zo>  x* ,  v* ,  O 
we  obtain  f*  by  minimizing 


G  m 


N 

7. 

i  =  l 


r  (zo,i>  '  8  >- 


(16) 


+See  equations  (2),  (4),  and  (6)  for  corresponding  boundary  relations. 
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CALCULATED 


3.81  CM.yCU.,  CONICAL  LINER.  CONE  ANGLE 


FIGURE  4.  PENETRATION  DEPTH  VS  VIRTUAL  STANDOFF  DISTANCE  FOR  THE 
3.81  CM..CU. .  CONICAL  LINER.  CONE  ANGLE  -  60*  . 


BRL  -  MR  -  1866 
CALCULATED 


FIGURE  5.  PENETRATION  DEPTH  VS  VIRTUAL  STANDOFF  DISTANCE  FOR  THE 
3.81  CM.,  CU.,  CONICAL  LINER.  CONE  ANGLE  «  90°  . 


20*  30*  40°  50°  60°  70°  80°  9 

CONE  ANGLE  ,  degrees 

FIGURE  6.  CALCULATED  VIRTUAL  ORIGIN  VS  CONE  ANGLE  , 


FIGURE  7,  CALCULATED  Umin  VS  CONE  ANGLE 


The  value  of  £  =  £*  which  minimizes  (16)  can  then  be  used  to 

determine  the  energy  constant  C  whenever  V.°  and  d.  are  known.  Since 
o  3  3 

both  dj  and  can  be  determined  with  BASC  then  a  ''best"  value,  C*, 

of  the  energy  constant  can  be  determined. 


VI.  SUMMARY 


We  have  shown  how  penetration  performance  -  stand-off  data  and 
hole  volume  -  stand-off  data  can  be  utilized  to  determine  values  of 
specific  functions  of  the  DiPersio,  Simon,  Merendino  shaped-charge 
parameters  C,  Umin,  and  t,,  and  that  specification  of  the  initial  jet 

tip  velocity  V.°  determines  "best"  values  of  t,  and  U  .  .  If,  in 

addition,  the  jet  diameter  d.  is  known,  then  the  energy  constant  C 

1  o 

is  readily  determined.  It  is  of  interest  to  note  that  V.  and  d.  are 

3  3 

readily  determined  from  Harrison's  BASC  code  and  are  functions  of  the 
liner  thickness  and  cone  angle.  The  implication  of  this  is  that  since 
x*  and  y*  are  determined  from  penetration  performance  data  one  may  then 
search  for  "best"  values  of  cone  angle,  a,  and  liner  thickness,  e, 
which  maximizes  the  jet  break-up  time  tj.  From  the  definition  of  x  we 
see  that  one  should  choose  a  and  e  such  that  V.  is  a  minimum. 
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